Ferromagnetism in the one-dimensional Hubbard model with orbital degeneracy: 

From low to high electron density 

Harumi SakamotoE, Tsutomu Momoi 
Institute of Physics, University of Tsukuba, Tsukuba 305-8571, Japan 

Kcnn Kubo 

Department of Physics, Aoyama Gakuin University, Chitosedai, Tokyo 151-8512, Japan 
(Dated: ) 

We studied ferromagnetism in the one-dimensional Hubbard model with doubly degenerate atomic 
orbitals by means of the density-matrix renormalization-group method and obtained the ground- 
state phase diagrams. It was found that ferromagnetism is stable from low to high (0 < n < 1.75) 
electron density when the interactions are sufficiently strong. Quasi-long-range order of triplet su- 
perconductivity coexists with the ferromagnetic order for a strong Hund coupling region, where 
the inter-orbital interaction U' — J is attractive. At quarter-filling (n = 1), the insulating ferro- 
magnetic state appears accompanying orbital quasi- long-range order. For low densities (n < 1), 
ferromagnetism occurs owing to the ferromagnetic exchange interaction caused by virtual hoppings 
of electrons, the same as in the quarter-filled system. This comes from separation of the charge and 
spin-orbital degrees of freedom in the strong coupling limit. This ferromagnetism is fragile against 
variation of band structure. For high densities (n > 1), the phase diagram of the ferromagnetic phase 
is similar to that obtained in infinite dimensions. In this case, the double exchange mechanism is 
operative to stabilize the ferromagnetic order and this long-range order is robust against variation 

of the band-dispersion. A partially polarized state appears in the density region 1.68 ~ n ~ 1.75 
and phase separation occurs for n just below the half-filling (n = 2). 

PACS numbers: 75.10.Lp, 75.10.-b, 74.20.-z 



I. INTRODUCTION 

The Hubbard model has been employed for a long time 
as a standard model for metallic ferromagnetism of itin- 
erant electronsoHu. It is, however, not easy to show 
that the Hubbard model simulates a metallic ferromag- 
net. Many investigations revealed that a simple single- 
band Hubbard model on a hyper-cubic lattice may not 
exhibit ferromagnetism. Some additional features in the 
model seem to be necessary in order to stabilize ferromag- 
netism. One of such features is a special lattice or band 
structure. Lieb first showecLthat a half-filled flat band 
induces a net magnetizationD. Then Mielke and Tasaki 
introduced other examples of flat-band ferromagnetstffl. 
Tasaki proved, furthermore, that a class of models with 
half-filled lowest bands with finite band dispersions ex- 
hibits ferromagnetism for sufficiently strong intra-atomic 
rcpulsionaJ. Miillcr-Hartmann found ferromagnetism in 
a one-dimensional system where the band dispersion has 
two minimau. It was also shown by perturbational as 
well as numerical methods that partially..!! 11 ed-bands in 
such systems can realize ferromagnetismallJO'O. 

Another important factor realizing ferroma gn ji ptism 
was discussed to be degeneracy of atomic orbitalst^tllEj. 
Actually magnetic carriers of typical metallic ferromag- 
nets originate from the degenerate 3<i-orbitals, though 
they are mixed with s- or p-orbitals. In the presence 
of degenerate orbitals, Hund coupling favors local triplet 
spin configurations of two electrons occupying different 
orbitals at the same site. Following mechanisms were 
proposed which lead to bulk ferromagnetism from intra- 



atomic Hund coupling. 

Zener proposed the so-called double exchange mecha- 
nism originally for doped manganites, which weeks, when 
the electron density per site is not an integerEJEJ, e.g., 
the density of electrons per site (n) is between I and 
I + 1. If electrons are strongly correlated, each lattice 
site is occupied by either / or / + 1 electrons, and Hund 
coupling forces their spins to be parallel in an atom. The 
magnitude of the formed atomic spin is 1/2 or (J + l)/2 
for I < N d - 1, and N d - 1/2 or N d - (I + l)/2 for 
I > N d , where N d denotes the degree of degeneracy. An 
electron can hop from a site occupied by / + 1 electrons to 
a site occupied by I electrons. Since the hopping matrix 
is diagonal with respect to the spin indices, an electron 
can hardly hop between sites whose spins are directed 
antiparallelly, while it can freely propagate among sites 
with parallel spins. As a result, ferromagnetic spin align- 
ment that lowers the kinetic energy is favored. 

Another mechanism leading to ferromagnetism is most 
effective when the electron density per site is an integer, 
n = I, where / is an integeMvhich satisfies 1 < I < N d — 1 
or N d + 1 < I < 2N d - IB. Each site is occupied by / 
electrons in the strong coupling regime and hoppings of 
electrons occur as virtual processes. The intermediate 
state of a second order process has the lowest energy if 
the spins are aligned parallelly. This leads to an effective 
ferromagnetic exchange interaction that induces orbital 
ordering at the same time. Van Vleck argued that the 
ferromagnetic exchange due to virtual hoppings-might be 
the origin of the metallic ferromagnetism in Nili3. 

Though the mechanisms favoring ferromagnetism can 
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be qualitatively understood as above, it is far from triv- 
ial whether ferromagnetic long-range order occurs in bulk 
systems. We need to examine whether the models with 
orbital degeneracy really simulate ferromagnets by reli- 
able methods. For this purpose, we study the simplest 
model with orbital degeneracy, the Hubbard model with 
doubly degenerate atomic orbitals. The model is de- 
scribed by the Hamiltonian 

H J2 ( C La C ima + H.C.) 
i<j m=A,B er=T,l 

i m=A,B i cr=T,| 

+{U' - J) 22 ^ n iA<jn lB< j 
- J ^Z^A]^ai c \bi^b] + H.c.) 

i 

- j 'J2( c Iai c Iai c *bi c iBI +H.c), (1) 

i 

where c\ ma (ci mCT ) creates (annihilates) an electron at the 
site i with spin a (=1, J.) in the orbital m (= A, £?). Here 
we assumed that only the hoppings between the same or- 
bitals exist, and the Coulomb interaction works only be- 
tween electrons on a same site. The interaction parame- 
ters U , V, J and J' represent the intra- and the inter- 
orbital repulsion, the exchange (Hund-rule) coupling, and 
the pair hopping, respectively. Particle-hole symmetry 
allows us to study only the case with < n < 2. 

Quite a few previous studies exist on the doubly degen- 
erate Hubbard model. One must be careful in comparing 
them, since different assumptions were employed for the 
interaction parameters. Many assumed that J' = 0. In 
fact, these parameters are dependent on each other if we 
calculate them from orbital wave functions. The relation 

J = J' (2) 

holds when the orbital wave functions are real. In addi- 
tion, the interaction parameters U, U' and J satisfy the 
relation 

U = U' + 2 J (3) 

for e g orbitals which is relevant for the cubic symmetry. 
This Hamiltonian has the spin SU(2) and charge U(l) 
symmetries. In addition, orbital degrees of freedom have 
U(l) rotational symmetry since 

[H,J2rf}=0. (4) 

i 

Here the orbital pseudo-spin operator is defined as rf = 
\ J2a=i i ^ia a i ^io-: where of denotes Pauli's matrix 
and Cj CT = {c iAa , c\ Ba ). The last symmetry holds if 
the coupling parameters satisfy the relations (^|) and 
(||), and it could be easily seen in the effective Hamil- 
tonian derived in the strong coupling limited. Because of 



this rotational symmetry, a rigorous argument for one- 
dimensional systems excludes-oscillating orbital LRO cor- 
responding to rf-correlationsE3. In the limit of J(= J') — 
0, the symmetry of orbital degrees of freedom becomes 
SU(2) invariant and the whole Hamiltonian has the high 
symmetry SU(4). 

We assume the relations (||) and (||) , and assume all in- 
teraction parameters to be positive throughout this work. 
In realistic situations the parameters satisfy U > U' > J, 
but we study the strong Hund coupling case J > U' 
as well, for theoretical interest. Let us explain the 
level structure of a doubly occupied atomic site since 
it rules the physics of the model ([[]). The spin-triplet 
states, where two electrons occupy different orbitals, have 
the lowest energy U' — J. In one of three spin-singlet 
states, electrons occupy different orbitals and the energy 
is U' + J. In the other two singlets with the energies 
U-J'(=U'+ J) and U + J'{=U' + 3 J), electrons oc- 
cupy the same orbital. In the strongly correlated regime 
where U' > J ^> almost all sites are empty or singly 
occupied for n < 1 , while they are singly or doubly occu- 
pied for n > 1 and spins in a doubly occupied site form 
a triplet. Thus the ground state properties for U' > J 
strongly depend on whether n is less than, equal to or 
more than unity. If J is larger than V, a doubly occupied 
site has negative interaction energy and hence electrons 
tend to form bound triplet pairs. In this case, sites are 
empty or doubly occupied for < n < 2 if N e is even. If 
N c is odd, there can be an unpaired electron, which will 
cause a strong finite-size effect. 

Roth examined the model in three dimensions in the 
quarter-filled case (n = 1)113. She found that the ferro- 
magnetic ground state has an orbital superlattice struc- 
ture in which two sublattices are occupied by electrons 
of different orbitals when the interaction is strong. This 
kind of ground state was studied by usin g . an , e ffec- 
tive Hamiltonian in the strong coupling limitta-EaEJ and 
by means of the mean field theory for general electron 
densityEi. Several exact diagonalization studies^eiie-DEL- 
formed for small clusters in one dimensionEjO'OmEj. 
Most of them investigated the quarter-filled case (n = 1) 

and found the ferromagnetic ground state for U' ~ J. 
Hirsch studied the less-than-quarter-filled case and found 
the ferromagnetic ground state. Their results, however, 
depend strongly on the boundary conditions and the 
number of lattice sites. There are rigorous proofs for 
the ferromagnetic groun d state in strong coupling limits 
in one dimensionr 3 rr 7 tr 8 l. One must be careful for that 
these proofs are valid in different limits of strong cou- 
pling. For the strong Hund coupling case (U — > oo and 
U' > J — > paL-pdstence of ferromagnetism is proved for 
1 < n < 2c3'EZl. In the special limit J = U' —t-po and 
U — > oo, ferromagnetism is stable for < n < 2cJ. Shen 
obtained a rather general result that the ground state is 
fully spin-polarized for any n between and 2 excep-t for 1 
if U = oo, and U'(> 0) and J = J'(> 0) are finiteM The 
present model (|l]) was also studied in infinite dimensions 
using the dynamical mean- field theorjCJCJ. Under the 
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assumptions (0) and (||), the ferromagnetic ground state 
with and without orbital order were found for n = 1 
and n > 1, respectively, while ferromagnetism was not 
found for n •f-ij— MP to the interactions of U = 40 and 
J < U' = 20eBEll. Biinemann et al. studied ferromag- 
netism in the three-dimensional two-band model with 
rather realistic DOS for Ni using a Gutzwiller approxima- 
tion. Their results showed that increase of DOS at the 
Fermi level could stabilize ferromagnetism for electron- 
doped cases (2 > n > 1)l3. 

In this paper, we studied the ground state of the 
model (0) in one dimension under open boundary con- 
dition, using the .-djensity-matrix renormalization group 
(DMRG) methodE-X This method enables us to study 
larger systems than those studied previously. Owing to 
the open boundary condition, we found remarkably lit- 
tle size-dependence. If we apply the periodic boundary 
condition, electrons can exchange their positions turning 
around the chain and it causes very large size-dependence 
of the ground state (e.g. even-odd oscillations). Prelim- 
inary results of this study were reported previously in 
Ref. H 

This paper is constructed as follows: In Sectio n p] , we 



briefly review the DMRG method. In Section [II, we 
study the present model with only nearest-neighbor hop- 
pings and report our numerical results, such as magnetic 
phase diagrams. Ferromagnetic phases appear for a wide 
parameter region and for all electron density (0 < n < 2). 
In most of density region, the ferromagnetic states are 
fully polarized, whereas partially polarized ferromagnetic 
states are found in a small region (1.68 < n < 1.75). 
We also found coexistence of metallic ferromagnetism 
and quasi-long-range order of triplet superconductivity 
when the Hund coupling is stronger than the inter-orbital 
Coulomb repulsion (J > U'). Near the half-filling case 
(1.75 < n < 2), phase separation occurs, and the param- 
agnetic (Haldane) phase and ferromagnetic one coexist. 
In Section IV, we examine the stability of the ferromag- 
netic states obtained in Section III, adding perturbation 



to the density of states. The metallic ferromagnetic order 
for more than quarter- filling (n > 1) is stable against this 
perturbation, but, for less than quarter-filling (n < 1), it 
easily becomes unstable by this variation. Finally, we 
give summary and discussion in Section |y|. Appendixes 
contain rigorous results in the strong coupling limit. Us- 
ing charge and spin-orbital separation, we show that the 
ferromagnetic spin state for less than quarter-filling is the 
same as one at quarter-filling. The ferromagnetism at less 
than quarter-filling is created by the exchange interaction 
produced by virtual second-order hopping processes. 



to White's paper for the DMRG method. We note only 
several points particular to this study. We employed the 
open boundary condition throughout this paper. In the 
DMRG calculations, we always kept rotational symmetry 
in the spin space so that the-pbtained ground state is an 
eigenstates of the total spintj. 

We used the finite-system method of DMRG to ob- 
tain the ground-state phase diagrams of the model with 
up to 16 sites. This method is superior to the infinite- 
system method in terms of accuracy of the ground- 
state quantities of finite systems. Another merit of 
this method is that the final result does not depend so 
strongly on the ground state of the initial small sys- 
tem and on the path of the renormalization process as 
in the infinite-system method. Repeated sweeps of the 
renormalization processes in the finite-system method 
can usually remedy the failure due to a wrong initial 
state. Near the phase boundary between a paramag- 
netic phase and a ferromagnetic one with saturated mag- 
netization, however, the final results sometimes depend 
on the choice of the renormalization paths. For ex- 
ample, two different choices of the paths (N, N c ) = 
(4, 2) (6, 3) (8, 4) (10, 5) (12, 6) (14, 7) (16, 8) and (N,N ) = 
(4,2)(6,4)(8,4)(10,6)(12,6)(14,8)(16,8) give a perfect 
ferromagnetic ground state and a paramagnetic one, re- 
spectively, at (N,N C ) = (16,8) for U' = J = 2|t|, where 
N c and N are the number of electrons and lattice sites, 
respectively. Since the transition is of first order, the 
energy levels of a paramagnetic state and a ferromag- 
netic one naturally exist very close to each other near the 
phase boundary. To study the phase diagram carefully, 
we searched the ground state through many renormaliza- 
tion paths. Finally, we regard the lowest energy state as 
the ground state. Hence, we need to try several paths 
and initial states and compare their results in order to 
obtain a reliable phase diagram as long as the number 
mo of reserved eigenstates of the density matrix is finite. 
The maximum number of mo was chosen to be 200 in 
this study for technical reasons. 

In order to calculate longer correlations we also studied 
the systems with up to 62 sites using the infinite-system 
DMRG method. Results of short chains were obtained 
by using the exact diagonalization method. 



III. RESULTS FOR THE NEAREST-NEIGHBOR 
MODEL 

In this section, we present the results for the model 
with only nearest neighbor hoppings (t{j = i#jj±i). 



II. METHOD 

We obtained the ground state of the model ([[]) with 
up to 62 sites using the-density matrix renormalization 
group (DMRG) methodEJ. The method is standard and 
we will not describe it in detail. Readers should refer 



A. Quarter-filled system 

We first exhibit the results in the quarter-filled case 
(n = 1). Figure [l] shows the (U'/t,J/t) phase diagram 
of the ground state of the quarter-filled system obtained 
for N — 4 and 8. The difference between the results for 
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N = 4 and TV = 8 is amazingly small. We hence be- 
lieve that the present results are representing the phase 
diagram of an infinite system, and did not study larger 
systems at this electron density. A ferromagnetic ground 
state with full spin polarization appears around J ~ U' 

for U'/t ~ 5. The parameter regions J/U' ^> 1 and 
J/U' <C 1 are paramagnetic. No partially polarized state 
was found at this filling. The present result is similar 
to those by previous studies with open boundary condi- 
tions though-different assumptions for parameters were 
employedESEj. In previous studies employing periodic or 
anti-periodic boundary conditions, the phase boundaries 
of their ferromagnetic phases depend on the system size 
greatljSlM 
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Here Si is a spin operator and Ti — (rf , rf , if ) is a 
pseudo-spin operator representing the orbital degrees of 
freedom, e.g., t 4 + is defined as ^ a cf Acr CiB<j- According 
to (^) the exchange energy between the sites i and i + 1 
is -2t 2 /(U' - J), -2t 2 /(U' + J) and -2Ut 2 /(U 2 - J 2 ) 



FIG. 1: Ground-state phase diagram at quarter-filling (n = 
1) of the model with nearest-neighbor hoppings. Interaction 
parameters are set as U = U' + 2 J and J = J'. The symbol 
"Stat = Max." denotes the perfect ferromagnetic phase, and 
"Stot = 0" the paramagnetic one. 



The appearance of the ferromagnetism in the region 
U' > J can be well understood in terms of the effective 
Hamiltonian in the strong coupling limit (U' — J ^> t). At 
quarter-filling, every site is singly occupied in this limit, 
and only the spin and orbital degrees of freedom remain. 
The second order effects of virtual hoppings between 
singly occupied sites lead to an effective HanMltanism Tt cS 
for the spin and orbital degrees of freedomliatjnJ'Ej: 



(5) 



(6) 



for the states |m, a)i\m, <r)i+i, \m, cr)i\fh, — o)i+i and 
|m, <j)i\m, — o)i+i, respectively, where \m, a)i denotes the 
state of the site i occupied by an electron of orbital m 
with spin a, and m labels the complement of m, e.g., 
A = B. As a result, H. c g favors ferromagnetism with an 
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orbital antiferromagnetic superlattice structure, where 
electrons on two sublattices occupy different orbitals. 
Numerical results of the four sites system described by 
H e ff exhibited that the ground state of 7i c ff is ferromag- 
netic for J > 0.35U'. This agrees with our observation 
that the lower boundary of the ferromagnetic phase in 
Fig. [I] approaches the line J = 0.35J7' for large U'/t. 
We may consider that the system is well described by 

HcH for U'/t ~ 10 and U' > J. A numerical study 
by two of present authors in infinite dimensions showed 
that the ferromagnetism with antiferromagnetic pseudo- 
spin i-order is stable for J > in the strong coupling 
limfO, which agrees with the result in the classical limit 
of Tleg- The ferromagnetic phase is smaller in one di- 
mension than that in infinite dimensions due to stronger 
quantum fluctuations. In one dimension, the antiferro- 
magnetic orbital correlations, which supports the ferro- 
magnetism, do not grow to a real long-range order but 
remains as a quasi- long-range order (QLRO), since or- 
bital degrees of freedom have U(l) rotational symmetry 
and a rigorous argument excludes oscillating orbital LRO 
in one dimensionLJ. 

In the ferromagnetic phase at n = 1, we found only 
fully polarized states. Note that the charge and or- 
bital degrees of freedom of fully polarized states can be 
described by a single-band Hubbard model with intra- 
atomic repulsion U' — J, where the conventional spin de- 
grees of freedom are replaced with the orbital oncsE3. 
Hence, we can conclude that this ferromagnetic state has 
orbital antiferromagnetic QLRO for U' — J > 0. In addi- 
tion, for general filling n, one can write down the asymp- 
totic form of the orbital correlation function in the fully 
polarized ferromagnetic phase in the form 



'o >i 



cos(tt.7tZ)|Z| 



(7) 



I 



for U' — J > 0. The exponent cts(= 1 + K p ) is obtained 
from the sipeie,-band Hubbard model with coupling U'—J 
at filling ttESEZI. We also checked the orbital ordering by 
the DMRG method and presented the results in Fig. |^. 

Also, from the above mapping, we can see that a metal- 
insulator transition occurs in the ferromagnetic phase on 
the line U' = J, i.e., the ferromagnetic state is insu- 
lating for U' > J and metallic for U' < J. (Figure 
shows chemical potential dependence of the electron 
density for U'/t = 15 and J/t = 10 with N = 16.) On 
the other hand, in the paramagnetic phase, since each 
band is quarter-filled and the perfect nesting no longer 
exists, a metal-insulator transition may not occur with 
an infinitesimal repulsion U' — J, which works between 
electrons with the same spin in different orbitals. Indeed, 
for paramagnetic states at J = 0, Assaraf et al. numeri- 
cally showed that a metal-insulator transition occurs at 
U c /t ~ 2.8iB. We expect that, in the paramagnetic phase 
with finite ,/(> 0), a metal-insulator transition may oc- 
cur at a finite positive value of U' — J. Figure || shows 
numerical results of the electron density as a function of 
the chemical potential in the paramagnetic phase with 
J/t = 3 and U'/t = 3, 4, and 8. The data for U'/t = 8 




FIG. 2: Correlation functions of orbital ordering ijiTj) in the 
ferromagnetic ground state for J/t = 15 and U'/t = 20 at 
filling n = 1 and n = 0.5. The system size is 34 and i is fixed 
at N/2 = 17. 



clearly show that the quarter-filled state is insulating. 
The data for U'/t = 3 and 4 seem to indicate that the 
systems are metallic at quarter filling, though it is hard 
to estimate the critical point from the present calculation 
because of finite-size effects. 

For J = 0, the effective Hamiltonian (|J) reduces to a 
model with SU(4) symmetry. It is known that the ground 
state of this model is a spin-singlet state (S = r = 0) 
in one dimension^. The paramagnetic state obtained 
for 0.3517' > J > in Fig. | may be interpreted as an 
extension of this SU (4) symmetric state. 

We find that, for the U' = J case, the ferromag- 
netic phase extends down to a weak coupling region. At 
V = J, double occupation of a site does not cost energy 
and the ferromagnetic ground state is metallic. The fer- 
romagnetism in the strong coupling region is stabilized 
by a kind of the double exchange .mechanism as was dis- 
cussed rigorously for U' = J = ocEl 

The ferromagnetic phase extends to the parameter re- 
gion J > U' . In this region, the fully polarized states 
can be-described by the attractive single-band Hubbard 
modeEa, whose ground state is known to have QLRO 
of the pairing correlationsc3. Hence the ferromagnetic 
ground state for J > U' has QLRO of triplet supercon- 
ductivity as well. For strong attraction (J —U' 3> t), all 
electrons are coupled into triplet pairs if N e is even and 
the system is described by an effective Hamiltonian 7ib 
in terms of hard-core bosons with spin unity: 
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defined as 
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FIG. 3: Electron density n per site as a function of the chem- 
ical potential fj, for U'/t = 15 and J/t = 10 for N = 16. 
The flat parts indicate insulating states. Abbreviations in- 
dicate the characters of the ground states as FM: ferromag- 
netic metal, FI: ferromagnetic insulator, FM(P): partially fer- 
romagnetic metal and PI: paramagnetic insulator. 



c 1 




FIG. 4: Electron density (n) per site as a function of the 
chemical potential (/x) for U'/t — 3, 4, and 8, and J/t = 3 for 
N — 16. The flat parts indicate insulating states. All states 
are paramagnetic. 
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where t = 2t 2 /{J - V) and J = 2t 2 /(J + U). The 
summation in 7ib runs over all the nearest neighbor sites. 
The creation and spin operators of a boson at a site i are 



L( t 



V2 



(9) 
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respectively. The number operator n\ is defined as 
S s =o ±ib\ s bis and P projects out states in which sites 
are occupied by more than one bosons. Since J is pos- 
itive, 7ib does not favor ferromagnetism and the ferro- 
magnetic phase should not extend to this region. We 
examined the phase boundary in the small system with 
N = 4 and found that the slope of the phase boundary 
approaches unity for large U'/t. The appearance of the 

ferromagnetism in the region J ~ U' is thus interpreted 
as a smearing out of the ferromagnetic ground state at 
J = U' to this region. Namely, this ferromagnetism is 
also caused by the double exchange mechanism at least in 
the strong coupling region. The above arguments should 
hold not only for the quarter-filled case but also for the 
hole-doped and electron-doped cases. Actually the slope 
of the phase boundary for J > U' approaches unity for 
large U'/t also in the systems with (N,N C ) = (6,4) and 
(4,6). Shen showed coexistence of ferromagnetism and 
triplet pairing correlations in. the whole region where 
U = oo and 1 < J/U' < ooEil On the other hand the 
present result shows that such a phase appears only in 
a small region close to the line J = U'. This discrep- 
ancy comes from the difference in the assumed relations 
between parameters. 

Next, we consider the spin-singlet ground state for 
J 3> U' where the system is described by 7ib- Figure 
H shows the density correlation function (n^n b j) of bosons 

and the spin correlation function (S^ ■ Sj) in the ground 
state of the system with 36 sites obtained by the infinite- 
system DMRG method. The correlations were measured 
from the center, i.e. i is fixed to be N/2. We observe an 
even-odd oscillation in the density correlations. This os- 
cillation is due to repulsive interactions between bosons. 
Actually the second term in 7ib works as repulsions since 
1 1 JS\ -S\ +1 + 2t-J\\ > holds when n\n\ +l = 1. On the 

other hand, {S\ ■ S^} shows an oscillation with the period 
four with slow decay. This oscillation is understood as a 
superposition of the density correlation with period two 
and the antiferromagnetic spin correlations caused by the 
exchange coupling J. 
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FIG. 5: Correlation functions (njn*) and (Sf -Sj) in the para- 
magnetic ground state for J/t = 50 and U'/t — 10 at quarter- 
filling. The system size is 36 and i is fixed at N/2 = 18. 

B. Hole- and electron-doped systems 

Next we study the hole-doped case (n < 1), where 
holes are doped into the quarter-filled system. Figure ^ 
shows the (U'/t, J/t) phase diagram of the ground states 
for n = 0.5 and 0.75 in the system with N = 16. A ferro- 
magnetic phase appears for a wide parameter region and 
the magnetization is saturated in the whole phase. Size- 
dependence among the results for N = 8, 12 and 16 was 
very weak, as in the case of quarter-filling. Comparing 
with the result for n = 1 we note that the ferromagnetic 
phase has expanded in the weak-coupling region down to 
U'/t ~ 3 as well as to a region where J > U' . 

We note that the lower boundary of the ferromagnetic 
phase for large U'/t (> J/t) almost agrees with that for 
quarter-filling. This result is a manifestation of the sep- 
aration of the charge and spin-orbital degrees of free- 
dom in one dimension. We can derive this result by us- 
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FIG. 6: Ground-state phase diagram for densities n = 0.5 and 
0.75 in the N = 16 system. Interaction parameters satisfy 
U = U' + 2J and J = J'. 



ing the same argument as that employed-by Ogata and 
Shiba for the single-band Hubbard modelcil. In the limit 
U' — J = oo j the motion of electrons is described by that 
of spinless fermions on a single chain, and the ground- 
state wave function is decomposed into a spinless fermion 
(charge) part and a spin-orbital part. The spin-orbital 
part is described in terms of spins and pseudo-spins in 
a squeezed system with iV c sites where empty sites are 
omitted. Interactions among the spins and pseudo-spins 
are caused by virtual hopping processes and described by 
the following Hamiltonian (see Appendix |Aj) 

, / sin2n7r\ 

ff = H 1 " -^r) Wcff ' (13) 

where Tt a s is the Hamiltonian given by Eq. (|^) where the 
site-labels of the spin-orbital quantities are those of the 
squeezed system. Since TC' e g differs from 7i s only by 
a multiplicative factor, it leads to the (U'/t, J/t) phase 
diagram identical to that for quarter-filling in the strong 
coupling limit. This argument explains the present result 
that the lower boundaries of the ferromagnetic phases for 

n = 0.5, 0.75, and 1 agree very well for U'/t ~ 10. We 
may conclude that the origin of the ferromagnetism in 
this region is the effective exchange coupling due to vir- 
tual hoppings. Orbital degrees of freedom are also simi- 
lar to those at quarter-filling except for the difference of 
the characteristic wavelength. This behavior can be seen 
in the orbital correlation function in the ferromagnetic 
phase obtained by the DMRG method (Fig. ||) and the 
asymptotic form given by Eq. (0). Note that the expo- 
nent as in Eq. (Q) does not depend on the filling n in the 
strong coupling limit. 

The above result greatly depends on the fact that the 
wave function is decoupled to the charge and spin-orbital 
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parts in the limit U' — J = oo, and that the ground- 
state energy is independent of the spin-orbital degrees 
of freedom in the first order of t. It should be noted 
that this is not the case in either higher dimensions or 
one-dimensional models with far-neighbor hoppings. In 
these cases the motion of electrons leads to effective spin- 
orbital interactions in the first order of i, as in the case of 
Nagaoka's ferromagnetism. In fact, ferromagnetism was 
not found for n < 1 in infinite dimensionscjO. We will 
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discuss the effects of far-neighbor hoppings in Section IV 



It is remarkable that the ferromagnetism appears for 
rather weak Hund-rule coupling, that is, J/t ~ 3 for 
U'/t ~ 5. Hirsch argued based on a diagonalization 
study of finite-size systems that the Hund-rule coupling 
is not effective enough to realize ferromagnetism in sys- 
tems with low density (n < 1) and that ferromagnetic 
exchange interaction between different sites are necessary 
to explain ferromagnetism in systems such as NicS. The 
present result, however, shows that a moderate Hund- 
rule coupling is sufficient to realize the ferromagnetism 
in the one-dimensional model with nearest-neighbor hop- 
pings. 

Finally we examined the electron-doped case where 
n > 1. The phase diagram for n — 1.25 in the systems 
with N = 8 and 16 is shown in Fig. 0. Size-dependence 
was very weak, as in the case of quarter-filling. In the 
systems with N c = Am + 2, where m is an integer, ground 
states in the weak-coupling region show small magnetiza- 
tion 5tot = 1 instead of Stot = 0. This weak spin polar- 
ization occurs when one-electron states at the Fermi level 
are half-filled, i.e. N c = 4m+2, under the open-boundary 
condition. The Hund coupling aligns the two electrons at 
the Fermi level parallel and produces the Stot = 1 ground 
state. We regard the appearance of Stot = 1 weak po- 
larization as a finite-size effect and not as an evidence of 
unsaturated ferromagnetism in bulk systems. In the fer- 
romagnetic phase, the magnetization is always saturated 
at this electron density. The ferromagnetism appears in 
a larger region than in the quarter-filled case in partic- 
ular for small J/t and large U'/t. This enhancement of 
the stable ferromagnetic state is a result of the double ex- 
change mechanism, which works effectively for 1 < n < 2. 
The appearance of ferromagnetism for U > U' > J ^> t 
is expected from a rigorous proof which holds in the limit 
where U > U' > J — > ooEZl. Our numerical results agree 
with this argument and assure that the ferromagnetism 
extends to a finite parameter region. 

We show the electron density n as a function of the 
chemical potential fj, for U'/t = 15 and J/t = 10 in 
Fig. H| for the system with 16 sites. For n > 1 and 
n < 1, the ferromagnetic ground state is metallic, since 
the compressibility dn/d/z is finite. For n = 1 and 2, 
the value of dn/dfj, is vanishing and hence the ground 
state is insulating. The system is a fully polarized fer- 
romagnetic insulator for the quarter-filling case. On the 
other hand, the ground state is paramagnetic at the half- 
filling (n = 2) where triplet electron pairs occupy all the 
lattice sites. Triplet pairs on nearest neighbor sites are 
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FIG. 7: Ground-state phase diagram at density n = 1.25 for 
the system with 8 and 16 sites. Interaction parameters satisfy 
U — U' + 2 J and J = J'. In the case of N = 8, paramagnetic 
states show small magnetization Stot = 1 because of a finite- 
size effect under open boundary condition. (See text.) 



interacting antiferromagnetically due to the second-order 
effect of virtual hoppings. As a result the ground state at 
n = 2 should be a Haldane singlet with an excitation gap. 
There is a jump in the density between n = 1.75 and 2 
in Fig. |[ which implies the occurrence of phase separa- 
tion. The occurrence r©f phase separation was found also 
in infinite dimensions^. 

We found partially polarized ground states in a small 
region of density where 1.6875 < n < 1.75. As the den- 
sity increases and approaches the half-filling, the anti- 
ferromagnetic interactions between the triplet pairs be- 
come operative and suppress the ferromagnetism due to 
the double exchange mechanism. This may cause the ap- 
pearance of the partially polarized states. We show the 
magnetization of a small system where (N, N e ) = (8, 14) 
as a function of J in Fig. |[ The figure shows that the 
partially polarized state appears in a finite range of the 
parameter, i.e., 8 < U'/t < 23 with J/t = U'/t-l. This 
result may imply that the transition from a paramagnetic 
to ferromagnetic ground state is of the second order for 

n ~ 2 in contrast to the discontinuous change to a fully 
polarized state observed in other lower density region. 



IV. EFFECT OF THE BAND-DISPERSION 

It is widely known that the profile of the density of 
states (DOS) of the band plays an important role in the 
realization of itinerant ferromagnetism. In the Hartree- 
Fock theory of the Hubbard model, for example, the fer- 
romagnetism is stabilized if U Dp > 1 where Dp denotes 
the DOS at the Fermi level. In fact, the whole profile of 
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FIG. 8: Total spin of the ground state as a function of U' /t 
in the system where J/t = U' /t - 1 and (N, N a ) = (8, 14). 



the DOS affects the stability of the ferromagnetisnxll 
The DOS of the nearest-neighbor model in one dimen- 
sion, which we have studied above, has a minimum at 
the center of the band and diverges at the band-edges. 
These features are quite different from those of DOS in 
higher dimensions, fn this section, we examine the ef- 
fect of the band-dispersion (or DOS) on ferromagnetism, 
studying a one-dimensional model with far-neighbor hop- 
pings. The model we employ has a linear dispersion of 
the band as 



-2t + (Atk/ir) (0 < k < tt) 
-2t- (itk/n) (-tt < k < 0) 



(14) 



instead of —2tcosk of the nearest-neighbor model, and 
hence its DOS is independent of the band energy (see Fig. 
|9|). For this model, the hopping integrals are calculated 
from the dispersion (14) by 



U 



N - 



— ^e fe sin(fci) sin(fcj), 

k 



(15) 



where the Fourier transform of c ima is given by 

N 



N - 



- sin(fci)4„ 



(16) 



i=i 



and k = £n/(N + 1) with I = 1,2,..,N. 
electron-hole symmetry of the dispersion dl^ 



Due to the 
) the system 



is bipartite, i.e. the hopping integrals tij vanish for \i — 
j\ = 2£ with integer £, and are always negative for = 
21+1 for positive t. The hopping integral t^ is given by 
-t[l-(-iy- j ]/ir(i-j) 2 in the limit N -> oo. We call this 
model linear-band model in the following. We study small 
clusters with N = 4, 5 and 6 under the open boundary 
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FIG. 9: Density of states per site and spin in the nearest- 
neighbor model (dotted line) and in the linear-band model 
(solid line). 



conditions, exactly diagonalizing them, and compare the 
results with those of the nearest-neighbor model. 

We show the ground state phase diagram for the 
quarter-filling obtained from the systems with (N, N e ) = 
(4, 4) in Fig. [Io| The phase boundaries for two mod- 
els are in good agreement in the strong coupling region 
where U' > J, while the ferromagnetism is slightly sup- 
pressed in the linear-band model in the weak-coupling 
region where J > V . The ferromagnetism is assisted 
by antiferromagnetic orbital QLRO when U' > J. The 
long-range hoppings in the linear-band model do not de- 
stroy the orbital correlations since Uj satisfies the bipar- 
tite condition. The present result suggests that the phase 
boundary between the ferromagnetic and paramagnetic 
phases does not depend strongly on the spatial depen- 
dence of the effective exchange interactions. 

Next we show the phase diagram for (N,N C ) — (6,4) 
in Fig. [ll] as an example of the hole-doped case. It is 
striking that the ferromagnetic region for the linear-band 
model is strongly suppressed compared to that for the 
nearest-neighbor model. The ferromagnetic phase ap- 
pears only for U'/t > 25 and exists in a very narrow 
region along the line J/U' = 1. This strong reduction 
of the ferromagnetic phase for U' > J is interpreted 
as a result of the breakdown of the separation of the 
charge and spin-orbital degrees of freedom. An electron 
can pass the other electrons by far- neighbor hoppings, 
and they lead to effective interactions among the spin 
and orbital degrees of freedom in the first-order of tij. 
The effective Hamiltonian in the limit U' > J — > oo, 
which contains only hopping terms and prohibits double 
occupancy, leads to a unique paramagnetic ground state 
(S = 0) in the linear-band model for (N,N C ) = (6,4). 
Consequently, the ferromagnetic phase does not appear 
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FIG. 10: Ground-state phase diagram at quarter-filling ob- 
tained by the exact diagonalization of the system with N = 
N c — 4. Filled circles and filled triangles display the phase 
boundary for the linear-band model and the nearest-neighbor 
model, respectively. 




FIG. 12: Ground-state phase diagram for (N,N C ) = (5,6). 
Filled circles and filled triangles display the phase boundary 
for the linear-band model and the nearest-neighbor model, 
respectively. 
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FIG. 11: Ground- state phase diagram for the system with 
(N,N e ) = (6,4). Filled circles and filled triangles display the 
phase boundary for the linear-band model and the nearest- 
neighbor model, respectively. 



for U' > J ^> t in this model and occurs only in the 
region J ~ U' and U'/t ~ 25. 

Figure [l2] displays the phase diagram for (JV, N c ) — 
(5, 6) as an example of the electron-doped case. In this 
case, the phase boundary for the linear-band model is 
similar to that for the nearest-neighbor model. The fcr- 
romagnetism for U' > J ^> t is caused by the double ex- 
change mechanism. Far-neighbor hoppings are expected 
to favor ferromagnetism as wel l as n earest-neighbor ones. 
We note that rigorous proofscjO'Es for the ferromag- 
netism in the strong coupling limit do not hold for the 
linear-band model. The present result implies that the 



ferromagnetic state for n > 1 is robust against the varia- 
tion of the DOS and would survive in higher dimensions. 
Actually, ferromagnetic state was found to be stable for a 
wide rangejXif interaction parameters for n > 1 in infinite 
dimensions^ . 

Though we have studied only small systems, the above 
results clearly show that the band-dispersion has a strong 
effect on the stability of ferromagnetism for the hole- 
doped system. On the other hand, it gives only a minor 
effect on the ferromagnetism for the quarter-filled and 
the electron-doped systems. 



V. SUMMARY AND DISCUSSION 

We studied the ferromagnetism in the one-dimensional 
Hubbard model with orbital degeneracy by numerical 
methods. We first studied the model with only nearest- 
neighbor hoppings by the DMRG calculations of clusters 
with up to 16 sites and obtained the following results: 

1. In the case of the quarter-filling (n = 1), a fully 
polarized ferromagnetic phase appears in a region 
U'/t > 4 and J/t > 4. The ferromagnetic state is 
insulating for U' > J and metallic for J > U' . The 
ferromagnetism for U' — J ^> t is caused mainly 
by the second-order hopping processes. The phase 
boundary between the ferromagnetic and paramag- 
netic phases approaches the asymptote J = 0.35(7' 
for large U' (> J). The ferromagnetic phase ex- 
tends considerably to a region with J > U' for in- 
termediate coupling regime, but the phase bound- 
ary on this side becomes parallel to the line J = U' 
for large U'/t. 



11 



2. In the hole-doped cases (n = 0.5 and 0.75), metal- 
lic ferromagnetism appears in a similar parameter 
region to that of insulating ferromagnetism at the 
quarter-filling. The phase boundary between the 
ferromagnetic and paramagnetic phases for large 
U'(> J) agrees with that for n — 1, which results 
from decoupling of the ground-state wave function 
into the charge and spin-orbital parts in the limit 
U' — J = oo. The ferromagnetic phase expands to 
the weak-coupling region more than the quarter- 
filled case. 

3. In the electron-doped case (n = 1.25), the ferro- 
magnetic phase dominates the region with U' > J. 
This ferromagnetism is metallic and caused mainly 
by the double exchange mechanism. Near the half- 
filled case (1.6875 < n < 1.75), this ferromagnetic 
order becomes partially saturated. 

4. In a large-Hund-coupling region, J > U', QLRO of 
triplet superconductivity coexists with the metal- 
lic ferromagnetism, which appears for all electron 
densities, i.e., < n < 1.75. 

It is important to note that the present results for the 
phase diagram showed little system-size dependence. We 
hence claim that the present results represent the prop- 
erties in the thermodynamic limit. 

The critical values of the interaction parameters for 
the appearance of ferromagnetism are shown in Fig. [l^. 
In the case of U' > J, the critical Hund coupling J c /t 
is almost constant for n < 1, while it suddenly drops for 
n = 1.25. Though we have not examined for other n 
values for n > 1, the data suggest a clear discrepancy 
between n > 1 and n < 1. For U' — J and U' < J, 
the critical inter-orbital repulsion U' c /t is apparently a 
smooth function of the density. For the case V = J, U'Jt 
seems to vanish for decreasing n. It would be interesting 
to examine whether the system exhibits ferromagnetism 
with infinitesimal interactions in the low-density limit. 

The double exchange mechanism causes ferromag- 
netism by lowering the kinetic energy, and hence the 
strength of the double exchange mechanism may be mea- 
sured by the difference AK of the kinetic energy in the 
ferromagnetic state Kp from that in the paramagnetic 
one Kp, i.e., AK — Kp — Kp. Largely negative AK im- 
plies that the double exchange mechanism is effective. 
We show AK for U' = J in the small systems with 
(N,N e ) = (6,4), (4,4) and (4,6) in Fig. [l|. We see that 
AK is a decreasing function of U'/t in all the cases and 
becomes negative for large U'/t. The value Ujt where 
AK = decreases with decreasing n which implies that 
the double exchange mechanism is more effective in the 
system with lower density for U' — J. The values of 
Ujt themselves are, however, much larger than the crit- 
ical value U'Jt for the appearance of the ferromagnetism. 
For example, U'Jt ~ 8 and U'Jt ~ 4 for (N, N c ) = (6, 4). 
This result shows that the ferromagnetism in the weak- 
coupling region is not caused by the double exchange 
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FIG. 13: Critical values of the interaction parameters for 
the appearance of ferromagnetism as functions of the den- 
sity: Jjt for U'/t = 20 (top), U'Jt for J/t = 20 (middle) and 
U'Jt (= J c /t) for U' = J (bottom). 
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FIG. 14: Difference AK of the kinetic energy in ferromagnetic 
states Kf from that in paramagnetic states Kp, i.e., AK = 
Kf — Kp . The numbers of system sizes and electrons are set 
as (N,NJ = (6,4), (4,4) and (4,6). 



mechanism. The ferromagnetism occurs for a gain of the 
interaction energy. The situation is similar in the system 
with (N, N c ) = (4, 6) for U' > J as well. The value J /t 
where AK = is much larger than the critical value J c /t. 
For example, J /t~7 and Jjt ~ 3 for U'/t = 20. The 
system-size dependence of UJt (and Jo/t) is left for a 
future study; whether it agrees with U'Jt (and Jjt) or 
not in the thermodynamic limit. 
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We found a partially polarized ferromagnetic state for 
a density region adjoining the phase-separation region 
just below half-filling, where the magnetization gradu- 
ally changes depending on the strength of interactions. 
In these states, ferromagnetic spin correlation decays es- 
pecially near edges of open systems. We may need fur- 
ther study in order to confirm that this partial ferromag- 
netism remains in the thermodynamic limit. There is also 
another possibility that a partially polarized state might 
apparently appear because of phase separation into a per- 
fectly polarized part and a paramagnetic part. Recently 
phase separation is found in other systems and.is.thought 
to be a common feature of strong correlationsHc3. 

Next, we examined the effect of the band-dispersion on 
the ferromagnetism using diagonalization of small clus- 
ters for a model with far- neighbor hoppings. The ferro- 
magnetism is greatly destabilized at less than quarter- 
filling n < 1. In this case, the second-order effects of 
hoppings, which are dominating in the nearest-neighbor 
hopping model, are overwhelmed by the first-order ef- 
fect. The phase diagram for the quarter-filling does not 
change greatly by far-neighbor hoppings that satisfy the 
bipartite condition. This result will be modified if the 
frustration is introduced by non-bipartite hopping inte- 
grals. The ferromagnetism for the electron-doped case is 
quite robust against variation of the DOS shape. 

We may extract some presumptions on the ferromag- 
netism due to the orbital degeneracy in two and three 
dimensions both frxan_the present results and those in 
infinite dimensionstSEil. At the quarter-filling (n = 1) 
the ferromagnetism is supported by the orbital antifer- 
romagnetic correlations for V > J. This mechanism is 
common in all dimensions and robust if the hopping in- 
tegrals satisfy the bipartite condition. We hence expect 
that the insulating ferromagnetic phase occurs in a simi- 
lar parameter region in all dimensions. The lower phase 
boundary will approach the line J = all' (0.35 > a > 0) 
for large U'/t, and a will decrease with increasing the 
dimensionality. For the electron-doped system (n > 1), 
it is natural to expect that the double-exchange mech- 
anism works in any dimension. The metallic ferromag- 
netism for U' > J will appear in a similar parameter 
region in any dimension for 2 > n > 1, since this fer- 
romagnetic order is insensitive to the variation of the 
band-dispersion. The situation is subtle in the case of the 
hole-doped system (n < 1). The double-exchange mech- 
anism hardly works for n < 1 because of inter-orbital 
Coulomb repulsion. In this density region, indeed, vari- 
ation of band-dispersion strongly destabilizes ferromag- 
netism in one dimension and no ferromagnetic phase was 
found in infinite dimensioris-jthough only sparse param- 
eter points were examinedoO. In three (and two) di- 
mensions, hence, the ferromagnetism might not appear 
in the simple nearest-neighbor hopping model on hyper- 
cubic lattices for n < 1. To confirm this argument, we 
need to take account of the interplay between the electron 
correlations and the band-dispersion and need to accom- 
plish serious calculations for two- and three-dimensional 



systems. 

For J > U', QLRO of triplet superconductivity ap- 
pears, which is supported by ferromagnetic order. In this 
region, a ferromagnetic state is stabilized by the double 
exchange mechanism. However, this ferromagnetic state 
is affected by far-neighbor hoppings and, in infinite di- 
mensions, the boundary of the ferromagnetic phase does 
not extend beyond the J — U' line up to the coupling 
U' = 1(M, though we (TM and KK) did not study the 
J > U' region. The coexistence of ferromagnetism and 
triplet superconductivity for J > U' might not appear 
in higher dimensions. We expect that a two-dimensional 
system may be a plausible candidate for showing ferro- 
magnetic triplet superconductivity. 

Properties of paramagnetic ground states are of theo- 
retical interest as well, though we did not study them in 
detail. The model has the S'/7(4) syjaametry for J = 
and solved exactly in one dimensioned. Recent numer- 
ical wo rks . clarified the properties of the quarter-filled 
systemMc3. It will be an interesting problem to investi- 
gate the spin and orbital correlations for general n and 
J > 0. From the charge and spin-orbital separation as 
shown in Appendix |X[ we can say that the spin degrees 
of freedom for n < 1 are equal to those at n — 1 in the 
strong coupling limit. In the other paramagnetic region 
for J ^> U', electrons are bound to form triplet pairs. We 
found strong charge correlations for the quarter-filling, 
but did not find LRO of charge density. It is an interest- 
ing problem to study whether a charge ordering occurs 
and coexists with antiferromagnetic LRO in higher di- 
mensions. 
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APPENDIX A: THE EFFECTIVE HAMILTONIAN 
FOR n < 1 IN THE STRONG-COUPLING LIMIT 

In this appendix, we derive the effective exchange 
Hamiltonian for n < 1 in the limit U' — J t in the 
nearest-neighbor hopping model. We treat the hopping 
term in ([!]) as a perturbation. In the unperturbed ground 
state, every site is singly occupied or empty. As the 
first-order effect of the hopping term, electrons hop from 
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singly occupied sites to empty ones. The effective Hamil- 
tonian in the first order of t is written as 



H. 



(1) 



err 



N 

;££(4 



(Al) 



where c ? - 



a 



.-.)(!• 



m — a ) IS 

the Gutzwiller-projected annihilation operator and hj is 
defined as ^ m a Cj mr7 Cj m<7 . In the unperturbed ground- 
state subspace hj is or 1. Since the electrons cannot 
exchange their positions under the open boundary condi- 
tions, the matrix elements of H^g are independent of the 

spin-orbital degrees of freedom. Hence is equivalent 
to a Hamiltonian for free spinless fermions, i.e., 



H. 



(i) 
err 



N 

^£(< 



H.c 



(A2) 



where a] is the creation operator of a spinless fermion. 
Then the number operator hi is equivalent to ajoj. We 
may write down the ground-state wave function | \T/ S ) as 
a direct product of the ground state of spinlessJjermions, 
|$sf), and a spin-orbital wave function |fl) aajj 



sf; 



(A3) 



Here |$sf) is a simple Slater determinant for spinless 
fermions and describes only the charge degrees of free- 
dom. The ground-state energy is degenerate for any spin- 
orbital wave function |0) up to the first order in terms of 
t. The second-order terms due to virtual hoppings deter- 
mine |f2). The effective Hamiltonian in the second order 
is given by 



H 
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err 



H 
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+ H 
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eff ' 



where 
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The first term H^ s °^ represents the exchange interac- 
tions between electrons occupying nearest-neighbor sites, 
where h(i, is defined in the effective Hamiltonian (|^) 
for the quarter-filling. The second term H^ffl represents 
correlated hoppings between next-nearest-neighbor sites. 

(2) 

The effective Hamiltonian H^ s ' contains both the charge 
and spin-orbital degrees of freedom. We may average 
charge degrees of freedom in H^J out since |$sf) is fixed 
to minimize HQ . In order to accomplish the averaging 
we introduce the squeezed system where the empty sites 
are omitted from the original system. Hence it is a chain 
with N c sites with spin and orbital degrees of freedom on 
each site. We label the sites in the squeezed system by rj. 
Employing the notation j v for the position in the original 
system of the 77-th electron in the squeezed system, we 
rewrite as 



and 



H 



[2b) 
eff 



N 

£££* 

i—l a rn 

1 est 



U' - J 



2+1 Trur 



H.c.) 



N. N 



H e« = E £ ^h^S^hiv, V + 1) (A7) 

r]— 1 i—l 



and 



N Q N 

^ff )= E££ £ H.c) a^efa, a , m ), (A8) 

77— 1 4—1 (Tq — ± mo — ± 
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where 



- + croS^ ) ( - - m T* J ( - + <?oS* +1 ) ( \ + to t* +1 



(t/') 2 - j 2 V 2 



j't 2 (\ \ (1 



C/ 2 - (J') 2 V 2 
J't 2 



" " CfCTo_roo Q—cra—ma 

JJ2 _ (Jl\2 1 V n+i »7+l 



U' — J \2 V V 2 

r/'/j 2 

(£/')2 _ J2 »7 '17 ^rf+l'ij+l 



{72 _ tji\2 \"n l 2 ? ' I n+1 I 2 ' 



(A9) 



1 

i-1 i i+1 



11 n+1 



FIG. 15: An examples of mapping from labelling of electrons 
by sites to labelling by counting only electrons from the left 



(77-th electron). A second-order term 77.jB4.cL 
mapped to the term t~ S^ +1 t+ +1 . 



The symbols <7o = +, — and too = +j — correspond to 
the spin j, | and the orbital A, B, respectively. The 
correlated hopping processes in the original system are 
mapped to the spin-orbital exchange processes in the 
squeezed system as illustrated in Fig. [l5|. 

In the limit where N — > 00 with n kept constant, 
the expectation values of 77,^+1(5^ and {nia^^i+i + 
R.c.)5ij 7i in the state |$sf) are independent of both 77 
and i. Then we obtain 



/ f7-( 2a )> 

( H cS , 



SF 



'hin i+ i)sF^2h(rj,r] 



Tj—l (To— ± mo— ± 



Here the notation (• • -)gF means the expectation value 
in the state l^gp). We note that the matrix elements of 
0(77, <T , to ) are equal to those of —h(r],r] + l)/2. We 
estimate the expectation values as 



\nin i+ i)sF = n 



(A12) 



and 



t \ sin2n7r / sinn7r\ ,..,„» 
{nial_ x a i+ i)$-F = n — — I — - — I , (A13) 

by using Wick's theorem. 

Finally, we obtain the effective Hamiltonian for the 
spin-orbital degrees of freedom as 



(-Hcff)sF = 77, 1 - 



sin 2?7,7r 
2nn 



Hefi- 



(A14) 



N, 



1) 



71=1 



Here Ti. c s is equal to the effective Hamiltonian (0) for 
the quarter-filling, where the summation over lattice sites 
are taken over those in the squeezed system. The factor 
77(1 — sin 277,77/277,77) is positive for all density (0 < 77 < 1). 
Hence, the effective spin-orbital exchange interactions in 
(AlO)the hole-doped system are same with that for the quarter- 
filled system except for a single multiplicative factor. The 
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FIG. 16: Ground-state phase diagram for (N,N C ) — (4,3) 
and (12, 9) as examples of the odd N e cases. 



present result suggests that the same fact holds quite 
generally for one-dimensional systems with any number 
of internal degrees of freedom. 



APPENDIX B: FERROMAGNETISM DUE TO A 
SINGLE UNPAIRED ELECTRON IN THE LIMIT 
WHERE J - U' = oo 

In finite-size calculations with odd number of electrons 
we find that the ferromagnetic state is stable in a large 
parameter region where J > U' as shown in Fig. |l6|. This 
ferromagnetism is caused by the motion of a single un- 
paired electron among bound triplet pairs. We consider 
a system with odd N c system in the strong coupling limit 
where J — > oo with J/U'(> 1) kept constant under the 
open boundary conditions. We prove, in the following, 
that the ground state of this system is fully polarized. 
Only nearest-neighbor hoppings are assumed to exist. 

In this limit two electrons on a site is bound to form 
a triplet pair. As a result M = (N e — l)/2 sites are 
doubly occupied and a single unpaired electron exits. In 
this limit, only allowed processes in the Hamiltonian are 
hoppings of an electron from a doubly occupied site to a 
singly occupied one and those from a singly occupied to 
a vacant one. The effective Hamiltonian H c n is written 
as 



N 
i=l 



T PfP i+1 C i+ i 



pD 



"Pi ' C im<jPi ' Pi+l c i+lmaPi+\ +H.C.}, (Bl) 



where Pf, Pf and P v are the projection operators for 
the doubly occupied triplet, the singly-occupied and va- 
cant states, respectively. We can divide the whole phase 
space into subspaces with fixed Na, Nb and S z . Here 
Na and Nb denote the number of electrons in the or- 
bital A and B, respectively (\Na — Nb\ = 1), which are 
conserved. Each subspace is further divided into n-\Cm 
sectors according to the configuration of doubly-occupied 
and vacant sites. States in each sector have the same con- 
figuration of doubly-occupied and empty sites if we ne- 
glect the singly-occupied site. Different sectors are dis- 
connected with respect to H e g. We call the sectors as 
Tj(Na, Nb, S z ), where j = 1, 2, ..,jv-i Cm- We can show 
for any two configurations a and j3 of spins and electrons 
belonging to a same sector Tj(Na, Nb, S z ) that: 

(I) there exists n > 1 such that ($ a \H%g:\$p) ^ 

(II) ($ Q |ff eff |^) <0. 

Here |$ a ) denotes a basis vector chosen as 



M 



ri(- 1 ^ 



A-OL \ ^Q2 ' ' ' A-q, n | 0) , 



(B2) 



where on denotes the atomic state of the z-th site and A ai 
is the corresponding operator, i.e. cIa-\ c Ibv c iAl c iB\> 



„t 



c\ai c \b-}) f° r doubly-occupied triplet 



states, c' ima . for singly-occupied states, and unity for 
a vacant state. The number ji denotes the position 
of the ^-th doubly-occupied site counted from an end. 
Here |0) is vacuum state. The proof of (I) and (II) are 
straightforwardEH. 

Now that the two conditions (I) and (II) on the ma- 
trix elements of H e g are satisfied, the Perron- Frobenius 
theorem ensures that the ground state in the sector 
Tj(Na, Nb, S z ) is unique and can be written as a linear 
combination of all the basis vectors with strictly positive 
(or negative) coefficients. The positive (or negative) def- 
initeness of the coefficients leads to that the ground state 
is an eigenstate of the maximum value of the total spinEH. 

We thus proved that the ground state in each sector is 
fully polarized. As a result the ground state for subspace 
with fixed Na and Nb is fully polarized. We understand 
that the large ferromagnetic region in the phase diagram 
in the case of odd N c for J — U' ^S> t is caused by the 
motion of a single unpaired electron and is a finite-size 
effect. 
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